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We analyse our recently-published velocity measurements of a Cen A (Butler 

et al. 2004). After adjusting the weights on a night-by-night basis in order to 

optimize the window function to minimize sidelobes, we extract 42 oscillation 

frequencies with I = 0-3 and measure the large and small frequency separations. 

We give fitted relations to these frequencies that can be compared with theoretical 

T^ , models and conclude that the observed scatter about these fits is due to the finite 

^ ! lifetimes of the oscillation modes. We estimate the mode lifetimes to be 1-2 d, 

O , 

Tj- , substantially shorter than in the Sun. 
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1. Introduction 

Its proximity and its similarity to the Sun have made a Cen A a prime target for 
asteroseismology. The first clear detection of p-mode oscillations was made by Bouchy & 
Carrier (2001, 2002) using velocity measurements with the CORALIE spectrograph over 13 
nights. These results confirmed the earlier but less secure detection by Schou & Buzasi 
(2000), who used photometry from the WIRE satellite. Based on the CORALIE data, 
Bouchy & Carrier (2002) published frequencies for 28 modes with / = 0-2 and n = 15-25 
and these have been compared with theoretical models by Thevenin et al. (2002) and Thoul 
et al. (2003). Meanwhile, oscillations in a Cen B have been reported by Carrier & Bourban 
(2003), adding to the growing list of stars in which solar-like oscillations have been detected 
- see Bedding & Kjeldsen (2003) for a recent review. 

We recently reported observations of a Cen A made over 5 nights in both Chile with 
UVES^ and in Australia with UCLES (Butler et al. 2004, hereafter Paper I). These obser- 
vations covered a shorter time span than those by Bouchy & Carrier but had substantially 
better velocity precision and the benefit of two-site coverage. 

In Paper I we described the method used to process the velocity measurements. In 
brief, we first removed slow drifts and sudden jumps of a few metres per second, presumably 
instrumental, from each time series. We then used the measurement uncertainties as weights 
in calculating the power spectrum (according to Wi = l/(jf), but we found it necessary to 
modify some of the weights to account for a small fraction of bad data points. In the end, 
we reached a noise floor of 2.0 cms~^ in the combined amplitude spectrum. In this paper, we 
extract oscillation frequencies, measure the large and small frequency separations and also 
estimate the mode lifetimes. 



2. Optimizing the window function 

The observing window contains gaps due to day time and bad weather, which means 
that each peak in the power spectrum is accompanied by sidelobes. This can be quantifled 
by calculating the spectral window, which is the power spectrum of the observing window 
and is shown in Fig. 6 of Paper I. The two-site coverage of our observations means that the 
gaps are relatively small in the latter part of the run, when the weather was best. However, 
as discussed in Paper I, the UVES data have much higher weight than the UCLES data. 



-'^ Based on observations collected at the European Southern Observatory, Paranal, Chile (ESO Programme 
67.D-0133) 
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due to their higher precision. This means that the spectral window, which takes the weights 
into account, has quite strong sidelobes (24% in power). As is well known, these sidelobes 
complicate the oscillation spectrum, especially for weaker modes, for which the interaction 
between noise and signal can cause sidelobes to become higher than the central peak and 
lead to mis-identifications. 

To address this problem, we have adjusted the weights on a night-by-night basis in order 
to optimize the window function. To be specific, we allocated separate adjustment factors to 
each of the three UVES nights and the five UCLES nights. The weights on each night were 
multiplied by these factors and the spectral window was calculated, and this process was 
iterated to minimize the height of the sidelobes. The final values for these weight multipliers 
were 0.88, 1.00 and 1.12 for the three UVES nights and 0.11, 0.42, 3.25, 8.19 and 7.62 for 
the five UCLES nights. We treated separately the short period of overlap (see Paper I), and 
the weight multiplier for that period turned out to be 1.39. 

This optimization process has significantly increased the weight given to the last three 
UCLES nights, which makes sense because these fill in the gaps in the UVES data. The 
power spectrum based on these adjustment factors is shown in Fig. 1. The resulting spec- 
tral window is shown in the inset, and the highest sidelobes are only 3.6% in power. The 
optimization process has slightly degraded the frequency resolution, with the FWHM of 
the spectral window increasing from 3.85 /xHz to 4.13 /iHz. However, the main trade-off is 
an increase in the noise floor in the amplitude spectrum at high frequencies, from 2.0 to 
2.9cms~^. Nevertheless, we have found that the improvement in the spectral window makes 
the oscillation spectrum much clearer, and this compensates for the reduction in signal-to- 
noise. We note that we found similar results to those reported here when we performed the 
analysis using a compromise weighting scheme, in which we used the square roots of the 
weight multipliers given above. 



3. Frequency analysis 

The large separation of a Cen A was measured to be ~ 106 /iHz by Schou & Buzasi 
(2000) and Bouchy & Carrier (2001). In Fig. 2 we show the central part of our power 
spectrum (2020-2970 /iHz) folded with a spacing of Az/ = 106.2 /xHz. We can clearly see 
peaks corresponding to modes with degree / = 0, 1,2 and 3. In Fig. 3 we show the power 
spectrum in echelle format, where we have smoothed in the vertical direction to make the 
ridges more visible. Again, four ridges corresponding to / = 0-3 are visible. 

The next step was to measure the frequencies of the strongest peaks in the power spec- 
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trum, which we did in the standard way using iterative sine-wave fitting. These frequencies 
are given in Table 1 and are shown as symbols in Fig. 4(a), over-plotted on a grey scale of the 
power spectrum itself (this time without any smoothing). The sizes of the symbols indicate 
the strengths of the peaks (the amplitudes are given in Table 2). We have assigned / values 
to the strongest peaks, based on their location in the echelle diagram. The four curves are 
fits to the data and are described below. 

There is a significant scatter of the individual peaks about each ridge, which we interpret 
as reflecting the finite lifetime of the oscillation modes. In several cases, the iterative sine- 
wave fitting has produced two peaks that appear to represent a single oscillation mode, which 
is to be expected if we have partially resolved the modes (see Sec. 4 for further discussion of 
the mode lifetimes). In these cases, both peaks are shown in Fig. 4(a) but in Tables 1 and 2 
(and also in Fig. 5) they are combined into a single weighted mean. It is also important to 
note that the separation between / = and / = 2 is not well resolved by our observations, 
and in one case (2572.7 /iHz) it seems that the two modes have combined to form a single 
peak. 

To show the power spectrum more clearly, we have plotted it on an expanded scale in 
Fig. 5 and indicated the extracted frequencies. Note that interference between noise peaks, 
oscillation peaks and their sidelobes has created apparent mismatches in a few places between 
the peaks in the power spectrum and the dotted lines representing the extracted frequencies. 
Such interference is relatively minor given the low level of both noise and sidelobes, but 
it is still significant (and it is, of course, the reason we chose to extract frequencies using 
the iterative fitting process). Also, as mentioned above, the iterative fitting has sometimes 
produced two peaks that have been averaged into a single value for this figure. 

Asteroseismology involves comparing observed frequencies with theoretical models, and 
the values listed in Table 1 could be used for that purpose. However, the scatter caused 
by the finite mode lifetimes means that we would do better to fit to the ridges, in order to 
obtain more accurate estimates for the eigenfrequencies of the star. We now describe the 
procedure that we used to perform this fit, which begins with an examination of the small 
frequency separations. 

Figure 6(a) shows our measured values for the three small separations, which are defined 
according to the usual conventions (for a recent review, see Bedding & Kjeldsen 2003). Thus, 
(5z/o2 is the separation between adjacent peaks with / = and / = 2 and 5vi^ is the separation 
between / = 1 and / = 3. The third small separation, (5z/oi, is the amount by which / = 1 
modes are offset from the midpoint between the / = modes on either side. Since one 
could equally well define 5voi to be the offset of / = modes from the midpoint between 
consecutive / = 1 modes, we have shown both versions in Fig. 6(a). 
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The observations show that the small separations in a Cen vary with frequency, as was 
already reported for 5z/o2 by Bouchy & Carrier (2002). We have fitted straight lines to each 
of the three separations (fitting a polynomial of higher degree is not justified by the data). 
These lines are shown in Fig. 6(a) and the equations are as follows: 

(1) 
(2) 
(3) 

Using these fits, we were then able to adjust the measured frequencies to remove the effect of 
the small separations and so align them all in a single ridge. The result is shown in Fig. 6(b). 
Note that we have taken the frequencies modulo Az//2, so that the even (/ = 0, 2) and odd 
(/ = 1,3) ridges are superimposed. To this single ridge we then fitted a parabola (the highest 
degree polynomial that is justified by the data), which is shown in the figure. 

This completes our fit to the observed frequencies, which has involved nine free param- 
eters (three linear fits and one parabola). We show this fit and the observed frequencies in 
Fig. 6(c), and also in Fig. 4(a). The equations for this fit are as follows: 

z/„,o = (2363.85 + 105.70 (n - 21) + 0.088 {n - 21)^) /iHz (4) 

z/„,i = (2414.35 + 106.10 (n - 21) + 0.088 (n - 21)=^) /iHz (5) 

z/„,2 = (2464.60 + 106.52 (n - 21) + 0.088 (n - 21)^) /iHz (6) 

Un,3 = (2509.70 + 106.41 (n - 21) + 0.088 (n - 21)^) /iHz. (7) 

These relations, rather than the actual frequencies, should be used when comparing with 
theoretical models because the measured frequencies are significantly affected by the mode 
lifetime. Note that these relations should not be extrapolated to frequencies beyond those 
actually measured. The formal uncertainties on the frequencies given by these relations are 
shown in Fig. 6(d), based on the uncertainties in the parameters of the parabola that was 
fitted in Fig. 6(b). Over most of the range the uncertainty is 0.3-0.4/iHz, rising to about 
0.8 /iHz at the ends. 

To examine whether the oscillations frequencies deviate systematically from equations (4)- 
(7), we have used those equations to "straighten" the ridges. Figure 4(b) shows the result, 
where we have summed the power of the four ridges (/ = 0-3) after first subtracting from 
each the frequency predicted by equations (4)-(7). Looking at the highest frequencies (above 
3000 /iHz), there seems to be evidence for a departure from the fit in the sense that the ridge 
of power curves to the left. More observations are needed to confirm this result. 

Another useful parameter to compare with models is the large separation, Au. The 
measured values are shown in Fig. 6(e) and there is an upward trend with frequency that 
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was previously noted by Bouchy & Carrier (2002) and which is responsible for the curvature 
in the echelle diagram. The line in Fig. 6(e) shows the large separation for / = 1 calculated 
from equation (5). Lines for other / values (not shown) are almost indistinguishable. 

The small solid circles in Fig. 6(c) show the 28 frequencies reported by Bouchy & Carrier 
(2002). Note these are all / = 0-2 because those observations were single-site and aliases 
prevented the detection of / = 3. On the other hand, our observations were shorter, which 
made it more difficult to separate / = from 1 = 2. Keeping in mind the scatter induced by 
the finite mode lifetimes, the agreement between the two results is excellent, especially for 
1 = 1. 



4. Mode lifetimes 

As already mentioned, the observed frequencies show a significant scatter about the 
fitted relations. If the modes were coherent (i.e., pure sinusoids with infinite lifetimes), 
we have shown using simulations that we would be able to measure them with a typical 
accuracy of 0.4 /iHz or better. The observed scatter is much greater than this, which we 
interpret as refiecting the finite lifetimes of the oscillation modes. As is well established for 
the Sun, the power spectrum of a stochastically excited oscillation that is observed for long 
enough will display a series of closely spaced peaks under a Lorentzian profile, the width of 
which indicates the mode lifetime (e.g., Toutain & Frohlich 1992). If, as in this case, the 
observations are not long enough to resolve the Lorentzian profile then the effect of the finite 
mode lifetime is to randomly shift each oscillation peak from its true position by a small 
amount. This is responsible for the scatter observed in Fig. 6. 

Measuring this scatter gives us an opportunity to estimate the mode lifetime in a Cen A. 
In Table 3 we give the scatter of the observed frequencies about the fit for two frequency 
ranges (a finer subdivision in frequency is not justified by the data). We have converted 
these to mode lifetimes using simulations, as follows. We simulated a single oscillation 
mode, sampled with the same window function and weights as the actual observations. We 
used the method described by Stello et al. (2004) to simulate an oscillation that was re- 
excited continuously with random kicks and damped on a timescale that was an adjustable 
parameter (the mode lifetime). For a given mode lifetime, we created 100 simulated time 
series that differed only in the random number seed. For each one we measured the position 
of the peak in the power spectrum and hence determined the scatter of the "observed" 
frequency about its true value. This was repeated for different values of the mode lifetime, 
allowing us to calibrate the relationship between the mode lifetime and the scatter in the 
observed frequencies. The inferred mode lifetimes for a Cen A are given in Table 3. 
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To check this technique, we also analysed segments of an 805-day series of full-disk 
velocity observations of the Sun taken by the GOLF instrument on the SOHO spacecraft 
(Ulrich et al. 2000). We sampled these data with the same window function and weights 
as the a Cen A observations, and extracted the frequencies of the strongest 30 oscillation 
peaks using the same iterative sine- wave fitting method as for a Cen A (see Section 3). We 
compared these frequencies with the values measured by Thiery et al. (2000) from the full 805- 
day GOLF series, and thereby determined the scatter on our "observed" solar frequencies. 
We then used the calibration derived from the simulations, as described in the previous 
paragraph, to infer mode lifetimes. The results for two frequency ranges are given in Table 3. 
In this case, we were able to compare these mode lifetimes with published values. For this we 
took values measured by Chaplin et al. (1997) from fitting Lorentzian profiles to a 32-month 
power spectrum of ground-based solar velocities. We see that our method reproduces the 
published lifetimes within the uncertainties, confirming the validity of the technique. 

Our main result is that mode lifetimes for a Cen A are 1-2 d, significantly shorter than 
the values of 3-4 d found in the Sun. Is this result consistent with the results for a Cen A 
presented by Bouchy & Carrier (2002)? Their frequencies were measured from a time series 
of considerably longer duration than ours (12.4 d versus 4.6 d), so we would expect them to 
have lower scatter. This is indeed the case - we measure a scatter of about 0.8 //Hz from 
their frequencies. This is much greater than would be expected for coherent oscillations 
(about 0. 1-0.2 yuHz), but is consistent with the mode lifetime that we measured from our 
observations. 

Note that the mode lifetime r is related to the FWHM F of the corresponding Lorentzian 
via r = l/(7rF). Thus, if a Cen A were observed for long enough, we would expect each 
mode to generate a Lorentzian profile in the power spectrum with width F = 2-3 /iHz. This 
is less than the FWHM of 4 /iHz of the spectral window of our observations, which explains 
why such a profile is not apparent in our power spectrum. 

We have assumed that the scatter of the observed frequencies about a smooth trend is 
due to the effects of the finite mode lifetime. However, we should also consider that the actual 
eigenfrequencies of the star may have small deviations from regularity. For the Sun, this is 
indeed the case, but at a low level: the frequencies reported by Thiery et al. (2000) deviate 
from a parabola in the echelle diagram by 0.2-0.3 /iHz. If a Cen A possesses irregularities of 
this magnitude, their contribution to the scatter that we measured is very small. Correcting 
for this would lead us to reduce the values of scatter in Table 3 by about 0.02 /xHz, and to 
increase the estimates of mode lifetimes by a few percent. We conclude that such frequency 
deviations have negligible effect on our measurements of the mode lifetime. 



5. Conclusions 

We have analysed the two-site observations of a Cen A that were described in Paper 1. 
Our main results are as follows: 

1. By adjusting the weights for each night, we optimized the window function and reduced 
the sidelobes in the spectral window to only 3.6% in power. This allowed us to measure 
the frequencies of the strongest peaks in the power spectrum using iterative sine-wave 
fitting, without any ambiguities from daily aliases. 

2. We clearly detected several modes with / = 3, the first time this has been achieved for 
solar-like oscillations. 

3. We made fits to the measured frequencies and obtained the relations in equations (4)- 
(7). Given the significant scatter in the frequencies, we suggest that comparison with 
theoretical models is best done via these fitted relations. 

4. Based on the scatter of the observed frequencies, we inferred mode lifetimes in a Cen A 
of 1-2 d. These shorter-than-expected lifetimes present a challenge to theoretical mod- 
els (Houdek et al. 1999), and also raise concerns for asteroseismology because they 
degrade the precision with which oscillation frequencies can be measured. It is clear 
that more stars need to be observed before we can determine how this important 
parameter varies in the H-R diagram. 
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Fig. 1. — Power spectrum of the combined a Cen A velocity time series, with weights for 
each night adjusted to optimize the window function. The inset shows the spectral window 
(power in arbitrary units), with the frequency scale expanded by a factor of 5. 
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Fig. 2. — The power spectrum of a Cen A folded at the large separation, showing peaks 
corresponding to / = 0-3. 
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Fig. 3. — The power spectrum of a Cen A in echelle format, and smoothed shghtly in the 
vertical direction. The column of numbers on the right shows the n value of each order. 



-12- 



3000 



I 
i 



>. 2500 



c 

3 
a' 

0) 



2000 



1500 




10 10 



Frequency modulo 106.2 yitHz 



Fig. 4. — (a) The greyscale shows the power spectrum in echelle format (without smoothing) 
and the symbols show the peaks extracted by iterative sine- wave fitting: / = (squares), 
/ = 1 (diamonds), / = 2 (triangles), / = 3 (crosses) and other peaks (plus signs). The sizes 
of the symbols are proportional to the amplitudes of the peaks. The curves show fits to the 
frequencies, as described in the text, and the numbers up the right side show the n value 
of each order, (b) The power summed over the four ridges (/ = 0-3) after first subtracting 
from each one the frequency predicted by equations (4)-(7). 
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Fig. 5. — Power spectrum of a Cen A with mode frequencies overlaid. Same as Fig. 1, but 
on an expanded scale and with the frequencies in Table 1 marked by dotted lines. 
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Fig. 6. — (a) The three small separations as a function of frequency. Symbols connected 
by dotted lines show the measured values: (5z/o2 (triangles), 61^13 (crosses) and Suqi (two sets 
of squares). The three solid lines are lines of best fit. (b) The frequencies after adjusting 
by subtracting off the small separations, to align them in in a single ridge. The curve is 
the parabola of best fit. (c) The measured frequencies (unadjusted) in echelle format. The 
curves show the four parabolas defined by equations (4)-(7), and the small solid circles show 
the frequencies reported by Bouchy & Carrier (2002). (d) Formal uncertainty in the fit as a 
function of frequency, (e) The large separation as a function of frequency for the four values 
of /, where the line shows the separation for / = 1 calculated from equation (5). In panels 
(b), (c) and (e) the symbols indicate / values, with the same meaning as in Fig. 4. 
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Table 1: Oscillation frequencies for a Gen A (yuHz) 
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Table 2: Oscillation amplitudes for a Cen A (cms 
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18 


13.3 


24.9 


18.7 


24.8 


19 


34.0 


18.6 


21.3 


15.2 


20 


20.7 


31.3 


31.8 


14.4 


21 


21.1 


31.6 


27.1 


19.6 


22 


16.7 


18.5 


40.3 


18.9 


23 


40.3 


26.2 


13.9 


11.2 


24 


27.8 


30.5 


18.0 




25 




14.4 






26 


21.4 


14.4 


16.8 




27 




7.5 







-16- 



Table 3: Mode lifetimes in a Cen A and the Sun 





Frequency 
(/.Hz) 


Scatter 
(/iHz) 


Inferred 
lifetime (d) 


Published 
hfetime (d) 


a Cen A 


1700-2400 
2400-3000 


1.50 ± 0.23 
1.60 ± 0.25 


1 4+0-5 
^•^-0.4 

1 O+0.5 
^^•-3-0.4 




Sun 


2400-3050 
3050-3700 


0.95 ± 0.17 
1.31 ± 0.23 


4-01?:? 

2 0+^-'' 

^•"-'-0.6 


3.17 ± 0.15 
1.59 ± 0.09 



